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ABSTRACT 

We present the discovery of the totally eclipsing long-period (P = 771.8 d) 
binary system WOCS 23009 in the old open cluster NGC 6819 that contains both 
an evolved star near central hydrogen exhaustion and a low-mass (0.45M Q ) star. 
This system was previously known to be a single-lined spectroscopic binary but 
the discovery of an eclipse near apastron using data from the Kepler space tele- 
scope makes it clear that the system has an inclination that is very close to 90°. 
Although the secondary star has not been identified in spectra, the mass of the 
primary star can be constrained using other eclipsing binaries in the cluster. The 
combination of total eclipses and a mass constraint for the primary star allows us 
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to determine a reliable mass for the secondary star and radii for both stars, and 
to constrain the cluster age. Unlike well-measured stars of similar mass in field 
binaries, the low-mass secondary is not significantly inflated in radius compared 
to model predictions. The primary star characteristics, in combination with clus- 
ter photometry and masses from other cluster binaries, indicates a best age of 
2.62 ± 0.25 Gyr, although stellar model physics may introduce systematic uncer- 
tainties at the ~ 10% level. We find preliminary evidence that the asteroseismic 
predictions for red giant masses in this cluster are systematically too high by as 
much as 8%. 

Subject headings: binaries: eclipsing — open clusters and associations: individual 
(NGC 6819) — stars: low-mass — stars: distances 



1. Introduction 



Long-period eclipsing binary stars are treasures for stellar astrophysicists because they 
offer the chance to examine the characteristics of (effectively) isolated stars with a high degree 



(Andersen 


1991; 


Torres et al. 


2010) 



more of the stars in the binary has evolved significantly off of the main sequence, the radius 
becomes an age-sensitive quantity. And if the binary is a member of a star cluster, this gives 
us a means of tigh tly constraining the age of the cluster and testin g stellar evolution theory 



quite strictly (e.g. iBrogaard et al.ll2012l 1201 ll : Meibom et all 120091 ) 



Although there have been a number of extensive variabi lity studies of NGC 6819 ( IKaluzny Shara 



19881 : IStreet et al.ll2002l . 12003k 120051 ; iTalamantes et al.ll2010l ). no eclipsing binary with period 
greater than about 15 days had been identified in the cluster. This is in part due to a low 
probability of occurrence, as well as a low probability of detection in studies of short dura- 
tion. Even though m any binary system s have been identified in a long-term radial velocity 
survey of the cluster (IHole et al.l 120091 ). there had not been a program to look for eclipses 
among the long-period spectroscopic binaries. The Kepler mission changed this by making 
continuous photometric monitoring possible, and thanks to a special effort by the science 
team to study clusters (an effort led by S. Meibom), a large portion of NGC 6819 (as well 
as the open cluster NGC 6791) is being observed. 



lr This is paper xx of the WIYN Open Cluster Study (WOCS). 
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Hole et al.1 (120091 ) originally discovered WOCS 230 09 (c^nm = 1 9 h 41 m lg;848, 5 : 



2000 



+40°07'27 / .'55; also known as KIC 5024447 and A851; lAunerl 1 19 74( 1 to be a single-lined 
spectroscopic binary and likely member of the cluster. We first detected the system as 
an eclipsing binary in quarter 4 of Kepler data, and flagged it as unusual because of the 
long duration of the eclipse (~ 1.8 d). Because the system's photometry places it close to 
the cluster's main sequence turnoff, the length of the eclipse does not stem solely from a 
long period, but is also because the radius of the primary star is larger than that of an 
unevolved main sequence star. A preliminary analysis indicated that the eclipse occurred 
near apastron, strongly constraining the inclination of the system to be very near 90°, or 
else the primary eclipse would not have been observed. Although we have only observed 
three eclipses of the system so far, the rare qualities of this binary make it possible to derive 
precision measurements of the individual stars already. We therefore present our initial 
analysis of the stellar characteristics in this paper. 



Observational Material and Data Reduction 



2.1. Spectroscopy 



Spe ctroscopic obse rvations were obtained as part of the WIYN Open Cluster Study 
(WOCS; iMathieul |2000| ). A detailed descr i ption of the acquisition and reduction of WOCS 
observations can be found in iGeller et al.l ( 120081 ). but a summary is presented here. WOCS 
spectroscopy uses the WIYN 3.5-m telescope^ on Kitt Peak along with the Hydra multi- 
object spectrograph (MOS), which is fiber-fed and capable of obtaining ~ 70 stellar spectra 
simultaneously (with ~ 10 fibers devoted to sky measurements). The observations presented 
here used the echelle grating providing a spectral resolution of ~ 15 km s _1 . The spectra are 
centered at 513 nm with a 25 nm range in order to cover a rich array of narrow Fe absorption 
lines near the Mg I b triplet. Spectroscopic observations of our target were completed using 
1 hour integrations per visit that were split into three 1200 s integrations to allow for the 
rejection of cosmic rays. 

Spectroscopic image processing was done within IRAF. After bias and sky subtraction, 
the extracted spectra were flat-field, throughput, and dispersion corrected. Calibration of 
the spectra utilized one 200 s flat field and two bracketing 300 s ThAr emission lamp spectra. 
Radial velocities were determined from a one-dimensional cross-correlation with an observed 



2 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatories 
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solar spectrum that was corrected to be at rest. The velocities were subsequently converted 
to the heliocentric frame and corrected for fiber offsets present in the Hydra MOS. The 
measurement precision is 0.4 km s -1 for o bservations o f sing le- lined systems (such as WOCS 
23009) down to a magnitude V ~ 16.5. iHole et al. J2009h presented a summary of NGC 
6819 radial velocity measurements, although the present study makes use of more recent 
spectroscopic observations by the same group. We used a total of 31 observations spanning 
from June 1998 to September 2010. The phased radial velocity measurements are plotted in 
Fig. [TJ showing the eccentricity of the orbit. 



2.2. Kepler Photometry 



Since quarter 1 as part of the Kepler Cluster Study (IMeibom et al.ll201ll ). the Kepler Sci- 
ence Team has recorded a series of 20 large image stamps (100x20 pixels) that cover the 
center of NGC 6819. Together they form a 200 x 200 pixel square field of view covering 
about 13'25 on a side. The image stamps were taken in long cadence mode, corresponding 
to exposure times of approximately 30 minutes. Although there are eight quarters of data 
on NGC 6819 that are publicly available as of this writing, we will discuss quarters 4 (19 
December 2009 - 19 March 2010) and 9 (19 March 2011 - 27 June 2011), which contained 
individual eclipses of WOCS 23009. 

We originally detected the primary eclipse in quarter 4 via simple aperture photometry 
conducted using modified versions of the PyK^f] software tools. For later analysis we made 
use of the publicly available light curve produced as part of the exoplanet detection program. 
For quarter 9 observations, a public light curve was not initially available, and so we deter- 
mined our own light curve using a pixel mask from quarter 5 (corresponding to the same 
observing season with the star on the same CCD). For our final results, we made use of light 



curves derived fro m the Pre-search Data Conditioning (PDC) pipeline ( IStumpe et al.l 12012 



Smith et al.l 120121 ). but we verified that these were good representations of the eclipses by 
comparing to our own reduction of the raw pixel data. We briefly describe that effort below. 

It is well known that the photometer o n the Kepler spacecraft introduces significant 
systematic trends into the photometric data (jKinemuchi et al.l 120121 ) . To begin to address 
this, we applied the kepcotrend algorithm using the first three cotrending basis vectors 
determined from analysis of light curves on the same CCD. After this procedure, there were 
still systematic trends visible outside of eclipses. For example, there appears to be spot 
activity with amplitude comparable to the secondary eclipse, but with period of around 5 



3 http : //keplergo . arc . nasa . gov/PyKE . shtml 
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d. We therefore opted to fit low-order polynomials to contiguous sequences of observations 
(generally separated from each other by momentum desaturation events on the spacecraft), 
deriving the normalization for the light curves from the fit. (This spot removal was done for 
PDC light curves as well.) 

Before finally normalizing the light curve, we also corrected the flux for crowding. Due 
to the large pixel size (4" width) and full-width-half-maximum (FWHM) of the point spread 
functions (95% encircled energy in a 4 pixel diameter aperture) of stars in the Kepler field, 
it is difficult to estimate the fluxes from Kepler photometry alone for stars in the crowded 
cluster fields. However, WOCS 23009 resides in a relatively sparse portion of the cluster, and 
crowding measures (ratio of target flux to total flux in the optimal photometry aperture) 
for quarters 4 and 9 (0.919 and 0.928, respectiv ely) from the Kep ler Data Search at the 



Mikulski Archive for Space Telescopes (MAST) ( IBrown et al.l 120111 ) support the idea that 
contamination due to other stars is small. We corrected the flux normalization for each 
quarter of data for this "third light" using the crowding measures above, so that the influence 
of light from other stars on the Kepler photometry should be negligible. We discuss nearby 
stellar sources more in §2.31 The corrected Kepler eclipse light curves are shown in Fig. [2} 
with the median out-of-eclipse magnitude (m KePtmed ) subtracted. 



We calculated the center of the Q4 primary eclipse using the method of lKwee fc van Woerden 



(11956( 1. and found barycentric Julian date (BJD) 2455267.6009±0.0007. The long duration of 
the eclipse (especially ingress and egress) allowed a more precise measurement of mid-eclipse 
than for most other eclipsing binaries in the Kepler field. Because the Q9 secondary eclipse 
occurs closer to periastron, the duration was considerably shorter (about 1.1 d), and there are 
fewer observations within the ingress and egress phases. We found BJD 2455695. 144±0. 014. 



2.3. Ground-based Photometry 



In order to more strongly constrain the binary period and obtain color information, 
we obtained VI images using the Mount Laguna Observatory (MLO) lm telescope on four 
successive nights (18-21 April 2012) to identify a second occurrence of the primary eclipse. 
The system was observed to be in eclipse on 20 and 21 April. 

The images were processed using standard procedures to apply overscan corrections from 
each image, to subtract a master bias frame, and to divide a ma ster flat fiel d frame. The 
photometry was derived using the image subtractio n software ISIS (lAlardll2000l ). The general 
procedure is described in iTalamantes et al.l (120101 ). as is some of the older photometric data 
that were used to help identify the reference out-of-eclipse level. Fig. [2] shows photometry 



-6 - 



from the four nights in April 2012. Although the ground-based images were not taken at 
the precise bottom of the primary eclipse, the depths of the eclipse in V and Iq agree well 
with the Kepler eclipse depth. 

In order to assess the potential contamination of the Kepler photometry from the light 
of nearby stars, we examined our deepest, best-seeing ground-based image in the Ic band. In 
the vicinity of WOCS 23009, stars had FWHMs of about 3.7 pix or about 1'.'5. Fig. M shows 
two much fainter stellar sources within 5" of WOCS 23009. We modeled the point-spread 
function of the image using DA0PH0T, and found that the two stars are about 3.9 and 5.5 Ic 
mag fainter respectively. We also subtracted off WOCS 23009 using ALLSTAR in order to look 
for closer resolvable stars, but found no additional sources. While these two nearby sources 
contaminate the light from the binary in the Kepler data, they should have negligible impact 
on the ground-based photometry using image subtraction techniques. 



3. Analysis 
3.1. Cluster Membership 



A first important question that must be asked is whether WOCS 23009 is a bona fide 
member of the cluster. It is not trivial to establish membership for a single star in a cluster, 
and arguments are generally framed according to probabilities. From an observational stand- 
point, the strongest membership arguments can be made when there are strong contrasts 
between the properties of the cluster stars and those of field stars in the same direction. 
Therefore, we look at several independent membership criteria. First, WOCS 23009 is 4'.33 
or ~1.8r n from the center of NGC 6819, where r c is a King model value of the core radius 
( IKalirai et al.l 120011 ) . At this radius, a simple ratio of the spatial density of probable clus- 
ter members versus the total spatial d ensity yields a 44% membership probability. This 
probability is estimated from Fig. 3 in IKalirai et al.l ( 120011 ). It should be noted that this 
membership probability is only a lower limit because no attempt was made to account for 
the differences in the cluster and field distributions in the color- magnitude diagram (CMD). 
The CMD position of the binary system mostly reflects the characteristics of the primary 
star (as shown earlier), and so the fact that the photometry places it among other single 
cluster members is additional circumstantial evidence of membership. 

The cluster membership can be estima ted more quantit atively from the radial velocity 
distribution for an unbiased sample of stars (IHole et al.ll2009l ). Our fit to the binary's system 
velocity 7 (see Tabled is very close to the cluster mean radial velocity (2.338±0.019 km s _1 ) 
and well within the velocity dispersion a = 1.009 ± 0.019 km s _1 . Based on the formulas 
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given in iHole et al.l . we calculate a 93% membership probability for WOCS 23009. Finally, 
we have an independent estimate of astrometric cluste r membe r ship. Star WOCS 23009 is 
too faint to be prese nt in the proper-mo tion study by ISanderd (119721 ) or the asteroseismic 
membership study of IStello et al.l (l2011bl ). We have calculated new proper motions around 
NGC 6819 down to V ~ 21 (Platais et al., in preparation). This study makes use of selected 
old photographic plates in combination with recent CCD observations at the 3.6 m Canada- 
France-Hawaii Telescope (CFHT). The accuracy of proper motions for well-measured stars 
near the center of the cluster is ~0.2 mas yr _1 . For star WOCS 23009, the uncertainty 
of proper motion is 0.10 and 0.18 mas yr _1 along the right ascension and declination axes, 
respectively. Images of this star are isolated and optimally exposed on all available plates 
and CCD mosaic frames. This gives us confidence that the star's astrometry is excellent. 
The calculated astrometric membership probability for WOCS 23009 is 99%. Summarizing, 
the body of evidence overwhelmingly supports cluster membership of this star, which should 
be considered a true cluster member. As such, the properties of the stars contain information 
on the properties of the cluster. 



3.2. External Constraints on the Binary 



Because WOCS 23009 is a single-lined spectroscopic binary, we needed to place an ex- 
ternal constraint on the mass of the primary star to be able to solve for the secondary star 
mass and both stellar radii. In this respect, the binary's membership in the cluster allows 
us to utilize external information that is not available for binary stars in the field. Models of 
the secondary eclipse indicate that the secondary star contributes about 0.19% of the total 
system light in the Kepler bandpass, so that out-of-eclipse system photometry can be taken 
as accurately representing the primary star. Therefore, a comparison to the photometry 
of other well-measured cluster stars can be used to constrain the primary star's charac- 



terist i cs. We assembled calibrated CCD photometry from available sourc es (IKalirai et al. 



200ll : iRosvick fc Vandenberd Il998l : IHole et all 120091 ; ISkrutskie et al.l 120061 ). and we deter- 
mined shifts from matching stars that would pla ce their photometry on t he same magnitude 
zero point as ours in B, V, and I g ba nds. (See iTalamantes et al.l |2010| for th e comparison 



withlRosvick fc V anden berglll998l . and I Jeffries et al. 



2001 and Hole et al. 



2009 



20121 for comparisons with IKalirai et al. 



) Magnitude offsets were never larger than 0.066 mag. The results 
are shown in Table [TJ For the filters where comparison is possible, the measurements are 
consistent within 0.03 mag. 



Though the stars in other NGC 6819 eclipsing binaries do not bracket the characteristics 
of the primary in WOCS 23009, they are close enough to allow us to make a modest extrap- 
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olation to estimate the mass of its primary star. IJeffries et al.l ( 120121 ) analyzed ground-based 



photo metry and sp ectroscopy for the eclipsing binary system WOCS 40007 (also known as 



A259; lAunerlll974l ). and Table [2] summarizes the results of that study, while Figure S] shows 
the deconvolved CMD positions of its stars. The secondary star in the WOCS 40007 system 
undergoes a total eclipse, so we can accurately measure the difference in V magnitude be- 
tween the components of the binary and use this as a consistency check on our calculations. 
The observed value (AV^ = 0.740) agrees with the theoretical predictions for the observed 
masses of the WOCS 40007 stars. 

By forcing isochrones to fit the photometry of the two WOCS 40007 stars at the ap- 
propriate masses, we can estimate the primary star's mass for WOCS 23009 using the V 
magnitude of the primary. We opted to rely more on the isochrone fits to the secondary star 
of WOCS 40007 because it is closer in mass (M s = 1.086 ± 0.018M Q ) to the Sun, so that 
its interior structure should be understood somewhat better. In addition, because WOCS 
40007 has total eclipses of the secondary star and there is a faint te rtiary on a long per iod 



orbit, the photometry of the secondary star is more tightly dictated ( IJeffries et al.ll2012l ). If 
the primary star of WOCS 40007 had been used to pin down the isochrones, the mass of 
the WOCS 23009 primary would be systematically larger. However, our quoted uncertainty 
on the WOCS 23009 primary mass encompasses fits that would match the WOCS 40007 
primary. (In other words, the photometry for both stars falls within the ranges predicted by 
the same isochrone for the measured masses and 1 sigma uncertainties.) 

With this algorithm, the theoretical estimates of the mass are relatively insensitive 
to the exact age of the isochrone used as long as the isochrones approximately match the 
morphology of the cluster turnoff. The evolution of the star's radius picks up dramatically 
for stars that have evolved past the reddest point of the kink in the isochrone (see §4.2[) . If 
the mass of any such stars are measured in the future, that will provide a strong upper limit 
to the mass of the primary of WOCS 23009. 

The reader should note that this is a differential estimate of the star mass, and so 
will avoid many of the potential errors involved in fits to the CMD (even when distance 
modulus and reddening are somehow constrained). Realistic systematic errors of 0.1 dex in 
the absolute metallicity scale only change the estimated mass by about 0.01M o , and potential 
helium abundance errors probably have a similarly small effect — to first order, composition 
errors should affect all of the cluster stars to approximately the same degree. Also, while we 
do base our mass determination on isochrones that fit the turnoff region of the CMD the best, 
we are not assuming a cluster age at this stage. Depending on the physics in the isochrones 
used (see Figs. @]and[5]), the age could be different, but the observational constraints from the 
CMD restrict the allowed values of the primary mass. The estimates produced by different 
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sets of isochrones agr ee well, with Dartm outh ( IDotter et a 



200 8]) isochrones pro ducing the 



lowest mass, Padova (IGirardi et al.ll2000f ) and Yonsei-Yale (IDemarque et al.ll20Q4f) isochr ones 
predicti ng a value about 0.005M ^ higher, and Victoria- Regina (IVandenBerg et al.ll2006l ) and 
BaSTI (IPietrinferni et al.l 120041 ) isochrones the highest by about O.O15M . We have chosen 
the Dartmouth value as our preferred value because of the overall fit of those isochrones to 
the turnoff region and the inclusion of significant physics updates (see §4.2p . 

The biggest source of uncertainty in the primary mass estimate comes from the masses 
of the WOCS 40007 stars because they allow for a small range of acceptable alternative 
fits. For a single set of isochrones, fits that return masses falling within the la uncertainties 
on both of the WOCS 40007 components produce a full range for acceptable fits of about 
0.046M Q . Therefore, we estimate Mi = 1.468 ± O.O3OM using the midpoint and somewhat 
conservative la error bars (slightly more than half of the full range). We emphasize that 
the precision and accuracy of the primary star mass estimate will improve as the data on 
WOCS 40007 improves and as additional binaries in the cluster are used to further constrain 
the analysis. 

Surface temperature is somewhat important for constraining the limb darkening coeffi- 
cients to be used in the modeling, but is even more important for the distance modulus deter- 
m ination in £14.21 Us i ng the system photometry and the empirical color-temperature relations 
of lCasagrande et al.l (120101 ). we can estimate the temperature of the primary from the various 
optical and infrared colors. We corrected for an average reddening E(B — V) — 0.12 ± 0.01 
Jeffries et al. 2012 ) after transforming this to other colors using the York Extinction Solvei0 



(jMcCallll2004l ). Because we did not calibrate Rc in our photometry, we cannot be certain 



it has the same zeropoint, and therefore, we give less weight to colors involving that filter. 
In addition, we give more weight to colors with a combination of small uncertainty in the 
photometry and small uncertainty on the reddening (when translated from uncertainty in 
B — V). Based on this, our best estimate of the primary star temperature is 7\ = 6320 ± 150 
K. 



3.3. Binary Star Modeling 



To simultaneously model the groun d-based radial velocities a nd photometry and Ke- 
pler photometry, we used the ELC code (jOrosz fc Hauschildtll2000l ). As discussed later, we 
fitted the data in two different ways. In one, we used PHOENIX model atmospheres to 
describe the limb darkening of the stars and the variation of emitted intensity with emergent 



4 http : //www2 . cadc-ccda.hia-iha.nrc-cnrc .gc . ca/ community/YorkExtinctionSolver/ 
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angle. In that case, we fitted the binary with a set of 10 parameters: orbital period P, time 
of conjunction (primary eclipse) t c , velocity semi-amplitude of the primary star K p , system 
radial velocity 7, eccentricity e, argument of periastron u, inclination i, ratio of the primary 
radius to average orbital separation Ri/a, ratio of radii R 2 /Ri, and temperature ratio T 2 /Ti. 
In the second set of trial runs, we fit for one additional parameter: one of the two coefficients 
in the quadratic limb darkening law for the primary star in the Kepler filter band. Because 
of the poor time resolution of the secondary eclipse and the faintness of the secondary star, 
the limb darkening coefficient for the secondary star is poorly constrained and was not fitted. 

In our models using the quadratic limb darkening law, one of the t wo coefficient s for the 
primary star was held fixed at the value given by ATLAS atmospheres (jClaretl |2000| ) for the 
estimated effective temperature and gravity of the star, while the other coefficient was fitted 
as part of the solution. By fitting for one of the coefficients, the effects of s ystematic errors in 



the o ther can be mitigated because the coefficients tend to be correlated (ISouthworth et al. 



20071 ). For comparison, we ide ntified the best-fit solution using PHOENIX model atmo- 
spheres (IHauschildt et al.lll997l ). Because the model atmospheres describe the variation of 
emitted intensity with emergent angle, there is no need to assume limb-darkening coeffi- 
cients. However, if there are differences between the vertical structure of the atmosphere 
in the models compared to the stars, this will lead to systematic error in the fitted stellar 
characteristics (and particularly the radii). Indeed, there are slight mismatches between the 
observations and best-fit model when using the atmospheres. 

The quality of the model fit was quan tified by an ov e rall y 2 , and the minimum value 



was s ought first using a genetic algorithm (iMetcarfd 1 19991 : ICharbonneaul Il995l ; lOrosz et al. 



20021 ) to e xplore a large swath of parameter space, followed by Markov chain Monte Carlo 
modeling (ITegmark et al.l 120041 ) to explore alternate models near the minimum and to esti- 
mate the uncertainties in the binary model parameters. The quoted parameter uncertainties 
are based on the range of values that produce a total x 2 within 1 of the minimum value, 
which approximates a la uncertainty (jAvnil 119761 ). The same procedure was used to deter- 
mine uncertainties in the radii of both stars and the secondary mass. (Uncertainties are not 
tabulated for the runs utilizing model atmospheres due to the computing time needed to 
explore the parameter space.) 

Because the uncertainties for the measurements are used in the calculation of x 2 , it is im- 
portant these uncertainties are as realistic as possible. We therefore forced the measurement 
uncertainties to be consistent with the observed scatter around the best fit model. If the 
measurement uncertainties calculated during the data reduction process are underestimated 
(as is often the case), this would inflate the total x 2 value and lead to an underestimation 
of the uncertainties on the binary star model parameters. To be conservative, after initial 
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optimization runs for the binary star model, the uncertainty estimates for the radial velocity 
and photometric observations were scaled upward to return a reduced \ 2 value of 1 for each 
type of measurement. (In other words, the reduced \ 2 value for radial velocities was used to 
scale the radial velocity uncertainties, and so on.) After this scaling (by the square root of 
the reduced x 2 )i we computed new models to determine the binary model parameters. 

Figure [6] shows the results of exploring the multi-dimensional parameter space during 
the binary star modeling for the fits with a limb-darkening law. The reader should note that 
the panels show the full range of the parameter space covered by models having \ 2 within 5 
of the best fit model. More than 10 6 models were run to survey the parameter space. Figs. 
[1] and [7] show comparisons of the radial velocities and Kepler light curves with the best fit 
model. Table [3] shows the parameters of the best fits using the limb-darkening law and model 
atmospheres. Most of the parameters are consistent within the uncertainties. Because the 
limb-darkening law fit matches the eclipse data more exactly, we will use the measurements 
from that run in the discussion below. 



4. Discussion 

Because we needed to apply an external constraint on the primary's mass to solve for the 
secondary star mass and both radii, our solution of the binary is not completely independent 
of stellar models, and will potentially be subject to systematic errors. Before discussing the 
results below, we make some estimates of how severe these errors might be. 

Given the measurable mass function (/ = Mf sin 3 i/(M x + M 2 f) for this single-lined 
spectroscopic binary, changes in the choice of primary mass will affect secondary mass deter- 
minations like dM 2 /dMi = 2q/(q + 3) where q = M2/M1 is the mass ratio. In other words, 
the systematic error in the secondary mass will be about 20% the size of the systematic 
error in the primary mass for this binary, assuming all other uncertainties are negligible. 
The masses of both stars influence radius measurements mostly through the total orbital 
velocity V = v± + v 2 because timing uncertainties are very small for the Kepler observa- 
tions of the long eclipses for this system. Even though we do not have measurements of the 
secondary star velocities or the semi-amplitude K 2 , the effect on the total orbital velocity 
goes like w ^ wr) « f^g> = j^dil + ||) = For the WOCS 23009 binary, 

this implies that fractional errors in radii will only be about 30% as large as the fractional 
systematic error in the primary mass. 



So even with a significant systematic error on the primary mass, a good constraint on 
the secondary mass and radii for both stars can still be obtained. If the primary mass has 
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a systematic error of about 0.045M Q , a systematic error of only 1% in the radii will be 
introduced. 

Because errors in the primary star mass and radius are correlated, the error ellipse in 
the M — R plane is tilted with respect to the axes (see Fig. [8]). This makes the error 
ellipse somewhat more parallel to the isochrones, potentially reducing the age uncertainty. 
However, after examining the best fit models (ones with \ 2 within 1 of the minimum value), 
we find the tilt is approximately 25°, which has a minimal effect. 

Although systematic errors in M 2 and R2 are both correlated with Mi (potentially 
making the error ellipse more parallel to the model main sequence, as in Fig. [9]), we find this 
effect to be even smaller, producing a tilt of only about 9°. The error bars plotted in Fig. [9] 
therefore give a very good approximation to the edges of the la error ellipse. As a result, it 
is unlikely that a systematic error in the primary mass is responsible for substantial error in 
the secondary radius. Regardless, surface gravity should be even less sensitive to potential 
errors because the correlations will mostly cancel. As illustrated in the lower right panel of 
Fig. El any correlation of log gi with the primary mass appears to be minimal, meaning that 
it is likely to be reliable. 

These exercises underline the fact that the primary star radius as well as the secondary 
star mass and radius are relatively insensitive to the primary star mass estimate for the same 
reasons — that this is a totally-eclipsing binary. 



4.1. The Low-Mass Secondary Star 



Only a handful of field stars are measured to have M < O.5M with small (< 3%) mass 
and radius uncertainties, and all of these are in eclipsing binaries. As Fig. [9] shows, the stars 
in binaries that have been studied generally show radii that are significantly la rger than mod- 
els p redict (see the figure caption for references on individual systems, and lLopez-Morales 
20071 for a review). The hypotheses that are currently being debated are that the larger radii 
stem from either a higher level of activity (with convection inhibited by strong magnetic fields 
or fast rotation) or a higher metallicity. Metallicity is difficult to measure directly for low- 
mass stars due to the complexity of their sp ectra, but indications a re that a metallicity tr end 
is only evident among single low- mass stars ( )L6pez-Moralesll2007l ) . ICoughlin et al.l (120111 ) re- 
cently described preliminary evidence from Kepler binaries that s pin-up due to a companion 
in short period binaries is responsible for inflated radii. However, iFeiden fc Chaboyerl (120121 ) 
have also emphasized that age is an additional uncertainty in comparisons of observations 
and models because even fewer systems provide age information. 
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The WOCS 23009 binary star system has made it possible to precisely measure the 
characteristics of a low-mass star in the cluster NGC 6819. This is the first such star measured 
in an open cluster, and as a result of its membership, we have excellent metallicity and 
age information from other cluster members. The most recent spectroscopic measurements 
indic ate a slightly super-so lar composition for NGC 6819 — in the most detailed study to 
date. lBragaglia et al.l (120011 ) find [Fe/H] = +0.09±0.03 from three red clump stars, where the 
quoted uncertainty is the error in the mean (systematic uncertainties are probably larger). 
Another unique aspect of this system is its large orbital period (a factor of 10 larger than 
any other known low- mass eclipsing binary), which allows the secondary star to evolve like 
it is in virtual isolation. 

As seen in Fig. EH our measurements indicate that the secondary star is one of the 
few well-meas ured low-mass stars with a radius that is consistent with model predictions. 
The models of iBaraffe et al.l (119981 ) pass slightly outside our conservative mass uncertainty, 
although the agreement would improve if those models used a chemical composition con- 
sistent with the super-solar metallicity measur ed for NGC 6819. Our measurements are 
co mpletely consistent with more recent models ( iDotter et al.ll2008l validated in calculations 
by lFeiden et al.l 1201 if ) that use an updated equation of state for the stellar gas. 



The closest observational comparison, CU Cnc A (IRibasI 120031 ). has approximately the 
same radius (R = 0.4317 ± 0.0052i? o ) and lower surface gravity (log g = 4.804 ± 0.011) for 
a slightly smal ler mass (M = 0.4333 ± 0.0017M Q ). If CU Cnc A is a member of the Castor 
moving group jRiba&]l2003f) . this would provide a metallicity ([Fe/H]= 0.0) and age estimate 
(320 ± 80 Myr; lRibasll2003h . Then CU Cnc A and WOCS 23009 B would provide a test of 
the idea that metallicity is the cause of larger-than-predicted radii. Before that is a strong 
test though, the mass uncertainty on WOCS 23009 B needs to be reduced, and CU Cn c A's 
membership in the moving group needs to be verified (see iFeiden fc Chaboyerl |2012| for a 
critique). 

The number of system s having low-mass stars a nd both reliable age and metallicity 
information is very limited. IFeiden fc Chaboyerl ( 120121 ) identify just 3 binaries with reliable 
information of this type, and all three have short periods that probably produce problems 
with stellar activity effects. Two (KOI-126 and CM Dra) contain stars near 0.2M Q , while 
the other (YY Gem) contains stars near O.6M . While the CM Dra and YY Gem systems 
have stars with larger-than-predicted radii, the two low-mass stars in KOI-126 have radii 
consistent with models for an age of 4 ± 1 Gyr, like that of the evolved third star in the 
system. 



WOCS 23009 B therefore gives us a very interesting opportunity: the opportunity to 
test models with a low-mass star that has well-determined and precise mass, radius, age, 
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and composition. The star has a radius and gravity that appear consistent with model 
predictions, and it is more compact than its closest observational comparison, CU Cnc 
A, at the 1.4a level. The long period of the binary assures us that tidal interactions are 
very unlikely to induce activity in the secondary star, unlike the situation for CU Cnc A 
(P = 2.771 d). In contrast, the low- mass eclipsing binaries with the pr evious records for 



the largest orbita l periods (about 4 1 d) are LSPM J1112+7626 A and B fllrwin et all 120111 ) 



and Kepler 16B (IDoyle et all 120111 ); both are found to be significantly larger than models 
predict, but do not have age or (in the case of LSPM Jl 12+7626) metallicity constraints 
that allow us to test alternate hypotheses. 

Low levels of intrinsic stellar activity could be partially responsible for the small radius 
of WOCS 23009 B, but the activity level of the star will be difficult to characterize because 
of the large difference in luminosities between the components of the binary. Although 
spot activity can affect radius measurements in eclipsing binaries, spot activity on the low- 
mass secondary in this system should be a minimal influence because the secondary star 
contributes such a small proportion of the system light. Based on the mass and likely 
age of the primary star, spots probably affect the photometry minimally. However, future 
eclipse observations with Kepler will allow us to monitor the system for unexpected changes 
accompanying the star's activity cycle. 



4.2. The Cluster Age and Distance Modulus 

We first examine the mass-radius diagram as this is the most direct way of comparing the 
observations to models. Radii of evolved stars can make good age indicators because they can 
be measured to high precision in eclipsing binary systems and they generally avoid sources of 
systematic errors (like distances and reddening) that can influence other age determinations. 
The biggest concern for age determination using the WOCS 23009 eclipsing binary is the 
reliability of the primary star mass estimate. As we emphasized in §3.3[ the primary star 
mass estimate is effectively constrained by a mass-luminosity relationship for main sequence 
(or near main sequence) stars. In contrast to the situation for the secondary star, there is 
a significant correlation between the primary star mass and radius. However, we evaluated 
this effect at the beginning of the Discussion section, and we use an appropriately tilted 
error ellipse in our analysis. 

To make the connection between stellar radius and age, theoretical models are necessary, 
and systematic differences in the input physics will affect the inferred age. While it is possible 
to find good fits to the primary star of WOCS 23009 when the isochrones are forced to match 
the components of WOCS 40007 in the CMD (as in Figs. H]and|5]), the ages implied by 
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these fits show significant differences, indicating that there are important disagreements in 
the physics inputs. Table H] summarizes some of the more important assumptions (from our 
point of view) in presently available model sets. In addition to those considerations, there has 
been an important ne arly 50% downward revision to the rate for the CNO cycle bottleneck 
reaction 14 N(p, 7) 15 ( jMarta et al.ll2008l ). Because the strength of the CNO cycle strongly 
affects the extent of convection in stars with small central convection zones, this probably 
has a significant effect on inferred ages for this cluster. None of the current publicly available 
isochrone sets includes th is rate revision with the exception of the Dar tmouth isoch rones, 
which use an earlier value (jlmbriani et al.ll2004i ) that is within 10% of the lMarta et al.l value. 
Due to the inclusion of this reaction rate as well as helium and heavy-element diffusion, we 
have used the Dartmouth isochrones as our primary theoretical comparison. Our future work 
will use additional stars in other eclipsing binaries in NGC 6819 and will employ isochrones 
with updated physics from various theoretical groups. 

Fig. [10] shows the comparison in the M — R plane. For the purposes of illustration, 
we have used identical ages for all of the theoretical isochrones, and when available, the 
same [Fe/H]. The precision due to uncertainty in the measurements of the primary star of 
WOCS 23009 is approximately 0.25 Gyr, but the systematic differences between isochrone 
sets (due to differences in coded physics) are similar to (if not larger than) this. Based 
on the Dartmouth set, the preferred age is about 2.65 Gyr. This is consistent to within 
the measurement uncertainties of the age implied by the low er-mass stars in th e binary 
WOCS 40007, and previous i sochrone fits to the CMD (2.5 Gvr. iKalirai et al.lboOll ; 2.4 Gyr, 
Rosvick fc Vandenbergi Il998l ). The Victoria- Regina, Padova, and BaSTI isochrones imply 
younger ages (with the Padova and BaSTI results lower by more than 0.25 Gyr or 10%). 

As we have emphasized, even though the primary star mass is somewhat uncertain, this 
affects the radius measurement in a minor way, and so it is still possible (in concert with 
a T e g estimate) to compute the luminosity of the star. Along with a bolometric correction, 
the absolute magnitude and apparen t distance modulus can be c alculated. We calculated 
the bolometric correction in V from IVandenBerg fc Cleml (120031 ). and used it along with 
the model radius and photometric temperature estimate from §3.31 to find My = 2.70 and 
(m — M)y = 12.39±0.10. The largest contributor to the uncertainty in this distance modulus 
by far is the temperature. 

Using the same method for calculating the distance modulus of WOCS 23009 A (includ- 
ing the employment of photometric T e g estimates), we determine (m — M)y = 12.36 ± 0.10 
and 12.41 ± 0.11 for the two components of WOCS 40007. The distance moduli from the 
three stars agree within the uncertainties, and the averaged distance modulus from the binary 
star components [(m — M)y = 12.39 ± 0.06] is slightly larger than previous determinations 
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(12.30 ± 0.10 from red clump stars, iJeffries et al.ll2012l ; 12.35 from fits to older isochrones, 
Rosvick fc Vandenberall998l ). but consistent within the errors. The use of a consistent tem- 
perature scale for all three stars will tend to minimize scatter in the distance modulus values, 
although systematic error in the distance moduli will still be present if there is a systematic 
error in the temperature scale. The most likely cause of such an error is the reddening, 
although metallicity error could also contribute. Because of their systematic nature, we add 
the contribution from the reddening ((Te(b-v) — 0.02 produces about 70 K uncertainty in 
temperature and 0.05 mag in distance modulus) and metallicity (cr[ Fe / H ] = 0.10 dex produces 
30 K uncertainty in temperature and 0.02 in distance modulus) in quadrature with the error 
in the mean for the three stars (0.06 mag) to get our final estimate (m — M)y = 12.39 ±0.08. 

We remind the reader that we chose the lowest (Dartmouth) estimate of the primary 
star mass to use in our binary star analysis. If the mass was systematically higher, this would 
result in a slightly larger radius (due to correlated errors), a larger calculated luminosity, 
and thus a larger distance modulus. So the mass estimate is unlikely to be responsible for 
our larger measurement. 

Any comparison in the CMD using just the distance modulus and reddening to shift 
the theoretical isochrones will immediately be subject to issues relating to color- T e g trans- 
formations, but we do a preliminary comparison here. Again, we use Dartmouth isochrones 
as our primary theoretic a l cho ice, here because they (and the Victoria-R egina set) employ 



the I VandenBerg fc Cleml (120031 ) transformations. IVandenBerg et a 



(2010) found that those 



relations show good agreement with the empirical relations from iCasagrande et al.l (120101 ) 
and with recent MARCS atmosphere models. In the (V, B — V) CMD of Fig. [TTJ the Dart- 
mouth models provide a good fit to the binary components and to single stars near the 
turnoff. The models appear to be slightly offset to the blue relative to the cluster stars, but 
this shift can be accounted for by the uncertainties in the reddening and distance, and in 
color-temperature transformations. The age indicated by the CMD position of the primary 
star of WOCS 23009 is consistent with the determination from the M — R plane. 

An examination of Figs. H] and [10] shows that the WOCS 23009 primary star is at an 
interesting point in its evolution, shortly before its radius begins to change rapidly. This 
evolutionary acceleration occurs shortly after the red kink in the CMD isochrones when 
convection in the core of the star starts to shut off shortly before central hydrogen exhaus- 
tion. Comparing mode ls of a given age with and without convective core overshoot (e.g., 
Pietrinferni et al.l 120041 ) . the radius of a star is unaffected by core overshooting until it has 
evolved past the red kink. Because the primary star has not reached the cluster's red kink 
at (V, B — V) ~ (14.8,0.65), the age determined from the stellar masses and radii are not 
sensitive to the amount of convective overshooting. The morphology of the turnoff and sub- 
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giant branch in the CMD are sensitive to overshooting, however. So if we use the binary star 
results in concert with the CMD, we can put new restrictions on the amount of overshoot. 
Larger amounts of convective overshooting bring more fuel to the fusion region of the core, 
and allows more massive stars to reach greater ages and higher luminosities before they turn 
onto the subgiant branch. When comparing to the Dartmouth models, the subgiant branch 
appears to be about 0.1 mag too faint, so this might be addressed with a larger amount 
of convective core overshoot. While the amount of overshoot in the Dartmouth models is 
similar to what is used in most other isochrone sets (see Table H] for a comparison of the 
overshoot in units of pressure scale height Hp), the lack of a match to the subgiant branch 
may be related to the lower (but more recent) 14 N(p, 7) 15 reaction rate used in the Dart- 
mouth models. Because this reaction sets the rate of energy generation in the CNO cycle 
and because the CNO cycle is responsible for convective cores in these main sequence stars, 
the revised reaction rate probably will affect stellar models in a manner similar to reduced 
convective overshooting. Further work combining analysis of the eclipsing binaries and CMD 
will be important in addressing this physics issue, which affects the masses of stars reaching 
the giant branch. 

The constraints that can be placed on the turnoff mass are also important to the astero 



seismic work on red giant stars in this cluster (IStello et al.ll2010l . l2011aU bl; iMiglio et al.ll2012 



mp 

m 



Corsaro et al.ll2012l ). Asteroseismology using the observed quantities Au (the frequency sep- 
aration of overtone modes) and z/ max (the frequency of maximum power) holds great promise 
for deriving the characteristics of isolated stars, but model-independent validat ion of the re- 



sultin g radii and masses is still underway (in eclipsing red giants, for example; iHekker et al. 
2010l ) . The so-called direct method uses Au and u m3iX (along with T e g) to compute mass 



and radius from scaling relations. Examination of the scaling relation linking Au and the 
mean density (Au/ Au & ~ \J pj p ) using stellar models shows that the relation probably 



holds to the level of a few percent for a wide range of evolutionary states (jStello et all 12009 
White et al.ll2011l ). The scaling relation involying u rasx is thought to relate to the atmospheric 



acoustic cutoff frequency (see iMosser et all |2010| for some empirical validation), which in 
turn relates to surf ace gravity and effective temperature (u max / u meLXyQ ~ (g/g©)/ ^/T cS /T eSjQ ; 
Brown et al.lll99ll ). When tested again st stellar models, this relation also appears to hold 
within a few percent (IStello et al.ll2009l ). When these two relations are combined, the radius 
and mass of a star can be calculated if T e g is known, although the asteroseismic observ- 
ables almost e ntirely cancel eac h other in the expression for mass due to the correlation 



between them (IStello et al.ll2009l ). "Grid-based" methods (finding best fits to Au, u m3iX , T e ff, 



and [Fe/H] values in a grid of models) avoid some of these problems, but it is still found 
that masses are less precisely determined than radii and subject to greater systematic errors 
(IGai et al.ll201ll ). While asteroseismic radii can be tested using techniques like interferom- 
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etry ( IHuber et all 120121 ). validation of asteroseismic masses is still at an early stage due to 
a lack of calibration stars with independent mass determinations. (See iKallinger et al.ll2010 
for an example of a calibration using a heterogeneous sample of red giants.) 

Because Kepler asteroseismic results have been used to derive masses for NGC 6819 
giants, a comparison with our eclipsing binary masses can be ma de. Two Ke p ler as teroseismic 
results on stellar masses in NGC 6819 have been reported: iBasu et all (120111) found an 
average giant mass of 1.68 ± O.O3M using a grid-based method, while iMiglio et al.l (120121 ) 
found a consistent value of 1.61 ± O.O3M purely using 4 different asteroseismic scaling 
relations involving various combinations of Au, z/ max , L, and T e g. For comparison, the best- 
fitting Dartmouth isochrone predicts that the mass at the bluest point on the isochrone (at 
the blue end of the subgiant branch) is about 1.52M and the mass of stars at the tip of 
the giant branch is about 1.56M . When any of the sets of isochrones are forced to fit the 
component s of WOCS 40007 (as in Figs. |4]and[5]), the red giant masses are lower than the 
Basu et al.l asteroseismic value by between about 0.07 (Yonsei-Yale, Padova) and O.13M 
(BaSTI). While the BaSTI result is probably affected by the lower amount of overshoot used 
in their models, there does not seem to be a set of models that simultaneously matches the 
masses and ra dii of the ec lipsing binary components, the morphology of the turnoff in the 
CMD, and the IBasu et all asteroseismic mass determination for red giants. 



Although the IMiglio et al 



result is consistent with some models (Padova, Victoria- 
Regina, and Yonsei-Yale) to within the uncertainties, it is still slightly discrepant with the 
preferre d Dartmouth models. In addition, three of the four asteroseismic scaling relations 
used by IMiglio et al.l involve luminosity, and so require a distance modulus. The distance 
modulus used by them is smaller than we find for the eclipsing binaries by about 0.18 mag. 
Adoption of the binary star distance modulus would result in more luminous stars and larger 
asteroseismic mass estimates, so this would worsen the d isagreement. In the same vein, the 
distance modulus for NGC 6819 derived by IBasu et al.l [(m — M) = 11.85 ± 0.05] is also 
smaller than implied by the binary stars. 

At this point, our results are indicative of a systematic error in the asteroseismic masses, 
although it must be remembered that this ultimately rests on the masses derived from the 
WOCS 40007 system and the radius of the primary star in WOCS 23009. Additional binary 
systems in the cluster are currently being analyzed and their measurements will improve our 
ability to predict turnoff and giant masses. The constraints from eclipsing binaries in this 
cluster and others in the Kepler field will help provide a more solid basis for asteroseismic 
determinations of mass if disagreements like these can be resolved. 
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5. Conclusions 

We present the discovery of total eclipses in a long period (771.8 d) eccentric binary 
system that is a member of the open cluster NGC 6819. The primary star is somewhat 
evolved and significantly larger than a main sequence star of the same mass, allowing us 
to derive an age of 2.62 ± 0.25 Gyr for the cluster using published models with the most 
up-to-date physics inputs. From the three cluster stars in eclipsing binaries that have been 
measured so far, we find an apparent distance modulus (m — M)y = 12.39 ± 0.08, where the 
uncertainty comes from the weighted mean and systematic uncertainties in reddening and 
metallicity. 

The measurement precision for the characteristics of this binary star (with the possible 
exception of the star masses) will improve with additional spectroscopic and photometric 
eclipse observations. Two of the largest uncertainties in the current analysis involve the 
stellar temperatures and the primary star mass. Better determination of the stellar surface 
temperatures will significantly improve the distance modulus determination, and a more 
precise determination of the reddening (via Stromgren photometry, for example) would allow 
a consistency check using colors. The primary mass will be more difficult to address, but 
new isochrone sets with improved physics inputs [such as the 14 N(p, 7) 15 reaction rate, 
diffusion, and convective core overshooting] and flexible chemical composition selection will 
help improve the estimate made here. However, it should be remembered that a solution 
for this binary system would not have been possible without constraints provided by other 
members of the cluster. Even with the relatively large uncertainty on the primary star mass 
(about 2%), we are able to reach a precision on radius measurements of better than 1% for 
both stars. In turn, this allows us to examine other stellar astrophysics questions. 

Even with the difficulties introduced by the large luminosity contrast between the two 
stars, the cluster age precision is determined to better than 10%. The characteristics of 
the stars we have studied so far are not sensitive to core convective overshooting, but the 
morphology of the turnoff in the color- magnitude diagram is, and this fact should make NGC 
6819 an excellent calibrator for the theoretical description of overshooting. 

The secondary star in the binary falls in the low-mass regime at about O.45M , but it 
does not show signs of radius inflation, as most other stars in the mass range do. Although 
high metallicity has been postulated as a cause of the larger radii, the secondary star contra- 
dicts this — NGC 6819's metallicity is super-solar ([Fe/H] = +0.09). While the secondary 
star in WOCS 23009 has a seemingly normal radius and it now holds the record for the 
low-mass star in the longest period eclipsing binary, two other relatively long-period systems 
(P « 41 d) contain stars that do show inflated radii. Episodic stellar activity may help 
explain these observations, and monitoring of the radii of stars in this and other long-period 
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binaries would help address whether this affects the measurements in an unexpected way. 
Further reduction of the mass uncertainty for the secondary will make this star a stronger 
test of low-mass stellar models. 

Finally, we emphasize that the mining of the ensemble of information available for this 
cluster has just begun. Although our results indicate that there is still some disagreement 
between masses derived from binaries and those derived from asteroseismology of giants 
in this cluster, the joint use of both techniques shows the potential for improving stellar 
modeling by pointing out where our understanding needs improvement. With additional 
constraints on the mass of the primary star, whether they come from asteroseismology or 
other eclipsing binaries, the age determination for NGC 6819 is likely to improve as the 
result of additional observations from the ground and with Kepler. 
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Fig. 1.— Phased WOCS radial velocities for WOCS 23009, along with the best fit limb- 
darkening law model (red line). The lower panel shows the observed minus computed values 
with error bars scaled to give a reduced \ 2 — 1 ( see §3.3p . 
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Fig. 2. — Top row: Kepler photometry of the primary and secondary eclipse of WOCS 23009. 
Bottom rows: Ground-based Vic photometry of the primary eclipse from April 2012. 



Fig. 3. — An examination of a 10" section of a ground-based Iq image to look for objects 
that could contaminate photometric measurements of WOCS 23009 (center). Detected stars 
are at 5 o'clock and 8 o'clock relative to the binary. Each pixel is approximately 0'.'4 on a 
side. N is up and E is to the left. 
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Fig. 4. — Color-magnitude diagram for NGC 6819 (data from iKalirai et al.l l200lh with 
probable single-star cluster members selected for V < 16.5 using iHole et al. j2009h . The 
system photometry for WOCS 23009 (dominated by the primary star) is shown with +, and 
the deconvolved photometry of the e clipsing stars in W OCS 40007 are shown with ellipses. 
Isochrones a re from the Dartmouth JPotter et al.ll2008[ [Fe/H]= +0.09; 2.5, 2.6, 2.7 Gyr), 
Yonsei-Yale dPemaraue et ali l2004[ [Fe/H]= +0.09; 2.4, 2.5, 2.6 Gyr), and Victoria-Regina 



( IVandenBerg et al.l 120061 : [Fe/H]= +0.136; ages 2.4, 2.5, 2.6 Gyr) groups. The theoretical 
predictions for stars with masses equal to the WOCS 40007 stars are shown with x, and the 
uncertainties in mass are delineated above and below those points. 
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Fig. 5.— Same as for Fig. H except for Padov a (IGirardi et al.ll200(tt [Fe/H]= +0.09; 2.2, 2.3, 
2.4 Gyr), and BaSTI jPietrinferni et al.lbooi [Fe/H]= +0.06; 2.1, 2.2, 2.3 Gyr) isochrones. 
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Fig. 6. — x 2 values for best-fit models using a limb darkening law as a function of binary 
model fitting parameters. 



31 




Fig. 7. — Top row: Comparison of Kepler photometry of the eclipses with model fits using 
a limb-darkening law {red line) and PHOENIX model atmospheres {blue line). Bottom row: 
Observed magnitudes minus model predictions for the limb- darkening solution with error 
bars scaled to give reduced x 2 — 1- 




Fig. 8. — Limb- darkened models with primary mass within la of the value extrapolated 
from WOCS 40007 components, and total \ 2 within 4 of the best fit model. 
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Fig. 9. — Mass-radius plot for low-mass members of eclipsing binary systems having mea- 
surement uncertainties of less than 3%. The secondary star in WOCS 23009 is shown with 
an open circle, and me asurements of other systems were take n from the following refer- 



ences: KOI126 B and C. ICarter et all (l201lh: CM Dra A a nd B. iMetcal 



Cnc A and B.lRibasI (J2003J); LSPM J1112+7626 A and B. Ilrwin et all fcOllh : Kepler 16B 



■e et all (11996( 1: CU 



Dovle et all (12011( 1: and NSVS01031772 A and B.lLopez-Morales et all (12006( 1. Model values 
are shown for |Fe /H|= and age 2.5 Gyr from Baraffe et al. ( 1998h . and [Fe/H]=+0.09 and 
age 2.6 Gyr from botter et all (bood ). 
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Fig. 10. — Mass-radius plot for members of eclipsing binary systems near the cluster turnoff. 
The error ellipse for the primar y star in WOCS 230 09 has the highest mass, and measure- 
ments of WOCS 40007 are from I Jeffries et all J2012h . Models have ages of 2.5, 2.75, and 3.0 
Gyr (from bottom to top). The Dartmouth, Yonsei-Yale, and Padova isochrones use [Fe/H] 
= +0.09, while the Victoria-Regina set uses the nearest tabulated [Fe/H]= +0.136, and the 
BaSTI set uses [Fe/H]= +0.06. 



-35 - 



C\2 




^ ID CO i> 



A 



Fig. 11. — Same as Fig. H] except that isochrones are shifted according to the preferred 
E(B — V) = 0.12 and (m — M)y = 12.39. Isochrones are from the Dartmouth group 
( IDotter et al.l l2008t [Fe/H]= +0.09; 2.5, 2.6, 2.7 Gyr), and ages were chosen to produce 
agreement with the bright component of WOCS 23009 (marked by +). 
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Table 1. Out-of-Eclipse Photometry for WOCS 23009 (A851) 



Filter 


T10/2MASS 


K01 


RV98 


H09 


B 


15.714 ±0.010 


15.716 


15.706 


15.710 


V 


15.090 ±0.010 


15.088 


15.104 


15.098 


Rc 






14.735 




Ic 


14.327 ±0.009 




14.356 


14.331 


J 


13.834 ±0.023 








H 


13.584 ±0.027 








K s 


13.556 ±0.045 









Note. 



T10: 


Talamantes e 


t al. 


( 


20inh. 


Skrutskie et al. 


(20061). K01: 


Kalirai et al. 



Hole et al. 



( 120091 ). 



Table 2. Summary Data for WOCS 4 0007 ( A259) Binary Components ({Jeffries et al 
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Primary 


Secondary 


M/Mq 


1.236 ±0.020 


1.086 ±0.018 


r/r q 


1.399 ±0.007 


1.098 ±0.004 


V 


16.138 


16.878 


(B-V) 


0.612 


0.715 


Teff (K) a 


6250 ± 150 


5870 ± 150 



"•Photometric tempe rature estimates 



(jCasagrande et al.ll2010l ). 



Table 3. Best-Fit Model Parameters for WOCS 23009 (A851) 



Parameter Limb Darkening Law Atmospheres 



Mi/Mo 


1.468 ±0.030 (constraint) 


7\ (K) 


6320 ± 150 (constraint) 


Xk,1 


0.3187 (fixed) 




Vk,i 


0.31 ±0.02 




X K ,2 


0.4013 (fixed) 




xv,i 


0.3419 (fixed) 






0.3507 (fixed) 




Xl,l 


0.1800 (fixed) 




yi,i 


0.3596 (fixed) 




/i 1 \ 

7 (km s ) 


2.22 ±0.02 


2.21 


P (d) 


77i.788±8:858 


771.790 


tc 


2455267.6000 ± 0.0004 


2455267.6000 


i (°) 


89.971« 


89.974 


e 


0.2550 ± 0.0004 


0.2551 


"(°) 


288.84 ±0.04 


288.84 


K P (km s" 1 ) 


6.96 ±0.13 


6.96 


Ri/a 


0.004860 ± 0.000016 


0.004899 


R2/CL 


0.00097651°™ 


0.0009634 


R 2 / R\ 


4 q77 +o.oo6 

' ' -0.009 


5.085 


(R 1 + R 2 )/a 


0.005836 ±0.000020 


0.005863 


T 2 /T 1 


0.572 ±0.003 


0.569 


M 2 /M Q 


n 447+0-010 
u -++' -0.011 


0.446 


Ri/Rq 


2.136 ±0.014 


2.151 


R 2 j Rq 


0.4292 ± 0.0033 


0.4229 


log gi (cgs) 


3.945 ± 0.004 


3.937 


log g 2 (cgs) 


4.823 ±0.009 


4.834 
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Table 4. Summary of Relevant Inputs for Model Isochrone Sets 



Isochrone 






Diffusion? 


Overshoot 




(Z/X) Q * 


AY/AZ 


Y 


Z 


^ov/ Hp 


BaSTI 


0.0280 


1.4 


N b 


N b 


0.14 


Dartmouth 


0.0267 


1.54 


Y 


Y 


0.20 


Padova 


0.0268 


2.25 


N 


N 


- 0.23 


Victoria- Regina 


0.0268 


2.2 


N 


N 


c 


Yale-Yonsei 


0.0253 


2.0 


Y 


N 


0.20 



a Initial values are tabulated if diffusion is incorporated. 

b Only the solar calibration for the BaSTI models incorporates 
diffusion of helium (Y) and heavy elements (Z). 

c The Victoria-Regina models use a different algorithm for con- 
vective core overshooting, so it is difficult to directly compare to 
the other sets. 



